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Outflow dynamics of dust-driven wind models and 
implications for cool envelopes of PNe 
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ABSTRACT 

The density profiles of cool envelopes of young Planetary Nebulae (PNe) are reminis- 
cent of the final AGB outflow history of the central star, so far as these have not yet 
been transformed by the hot wind and radiation of the central star. Obviously, the 
evolution of the mass loss rate of that dust-driven, cool wind of the former giant in 
its final AGB stages must have shaped these envelopes to some extent. Less clear is 
the impact of changes in the outflow velocity. Certainly, larger and fast changes would 
lead to significant complications in the reconstruction of the mass-loss history from a 
cool envelope's density profile. 

Here, we analyse the outflow velocity v cxp in a consistent set of over 50 carbon-rich, 
dust-driven and well "saturated" wind models, and how it depends on basic stellar 
parameters. We find a relation of the kind of v exp oc (L/M) . By contrast to the 
vast changes of the mass-loss rate in the final outflow phase, this relation suggest only 
very modest variations in the wind velocity, even during a thermal pulse. Hence, we 
conclude that the density profiles of cool envelopes around young PNe should indeed 
compare relatively well with their recent mass-loss history, when diluted plainly by 
the equation of continuity. 

Key words: stars: AGB and post-AGB, stars: mass-loss, stars: winds, outflows, stars: 
supergiants, circumstellar matter, planetary nebulae: general 



> . 1 INTRODUCTION 



3 



The intriguing detail displayed by PNe in visual wave- 
lengths is the result of complex and dynamic interaction of 
very different phases of circumstellar matter, following rapid 
changes in stellar evolution. The ionizing UV radiation and 
hot, fast wind of the hot central star modifies a c ool outer 
envelope from inside-out (|Schonberner et al]|2005t ). In fact, 
many circular and elliptical PNe often show a double-shell 
structure which consists of a rim - seen as a bright inner 
ring - and the surrounding shell (jPhillips et al.ll2009l . and 
references therein). The latter is material already heated up 
by the ionizing radiation. The inner rim, on the other hand, 
which appears to move with a velocity of ~40 km s , is 
caused by the interaction of the hot, fast wind of the cen- 
tral star with the cool, slowly expanding shell. The latter is 
the product of a recently ceased, cool, slow (~15 km s _1 ), 
and dense dust-driven outflow, which accompanied the final 
asymptotic giant branch (AGB) stages of the star. 

Renzini once coined the nickname "superwind" 
ijRenzinil fl98lh for this final, massive and dust-rich AGB 
mass loss. It is a result of the very efficient interaction of the 
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stellar radiation of a cool AGB giant and the newly formed 
dust particles in its atmosphere. This superwind removes of 
the order of one solar mass in the course of about 30 000 
years from the stellar envelope, prior to the very end of the 
AGB phase. This concerns stars with low- and intermediate 
initial mass, of about (~l-8 Mq), and their mass-loss rates 
can reach from 10~ s Mq yr" 1 to 10" 4 Mq yr" 1 . 

A physical record of such a cool outflow can be found 
outside the visible structures of a young PN. By contrast to 
its spectacular optical counterpart, the cool envelope is not 
at all so easy to observe. Only the remarkable advances of 
observational means in the past two decades in the infrared 
(ISO satellite in the 90ies and Spitzer in the past decade) and 
millimeter-waves (IRAM, large advances in detector technol- 
ogy, lead ing to a p romising potential of upcoming ALMA, 
see e.g. |Coxl 1200 ll ) have put such cool envelopes now in 
reach of human scrutiny. 

Still, details of that outflow phase, and how it gives 
way to the observed variety of PNe - like circular, ellipti- 
cal, or bipolar - are not yet fully understood. It seems clear 
only that the transition occurs when the stars are develop- 
ing from the AGB to PNe, i.e. when they are in the short 
proto-planetary nebula phase l|Fong et alj I2OO6L Unfortu- 
nately, this phase is rare to observe as it is short-lived (of 
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the order of 10 3 years). But the younger a PN is, the less 
has the evidence of its recent past - imprinted in its cool 
outer envelope - been modified by the hot central star. 

In addition to all these questions of PN-formation, there 
is a further motivation to fully and quantitatively under- 
stand the final phase of mass loss on the AGB: The mass 
lost during this time determines, by a large proportion, how 
much mass actually remains in the stellar remnant, the white 
dwarf. In other words, work on superwinds als o aims to re- 
prod uce the observed initial-final mass relation l|Weidemannl 
1200(f ). 

For a quantitative understanding of the radial density 
profile of a young PN's cool, outer envelope, the mass-loss 
history is clearly not the only factor of importance. Equally 
well we should know, how fast the outflow was and how 
much its velocity was varying. Nevertheless, some cool PN 
envelopes appear to be reproduced quite well by means of a 
proper mass-loss history alone, by making the simplification 
that on the star's evolutionary time-scale the outflow veloc- 
ity did neithe r change too much n or too fast. This approach 
was taken bv lPhillips et al.1 (|2009t ). who compared observed 
density profiles of PNe derived from Spitzer images to theo- 
retical predictions of stellar evolution calculations which in- 
cluded a parameter-dependent mass-loss rate. The latter is 
based on models for a C-rich, dust-driven, and dense wind. 
Despite its simplicity, the assumption of constant outflow 
velocity still leads to reasonable agreement between models 
and observations. Hence, here we attempt a more detailed 
investigation of how much the outflow velocity is actually 
changing in the dense, C-rich, and dusty envelope of a red 
giant, close to the very end of its AGB stage. 

There are, of course, immense limitatio ns to the ob- 
servat ional evidence of such objects. Hence, lHabing et al.l 
l| 19941 ) attempted such an investigation of outflow velocities 
based on theoretical arguments. For high mass-loss rates, 
M > 10 -5 MQyr -1 , they found a very modest dependence 
on stellar quantities: v out oc L°' 3 (5 ' 5 (M/a) ' 04 , where L* is 
the stellar luminosity, S the dust-to-gas ratio, and a the grain 
size. The latter appears to have a nearly negligible influence 
on the outflow velocity. The same study found, observational 
evidence would match well such a small dependence of the 
outflow velocity on stellar parameters, i.e. on the luminosity. 

We here want to take a new and different look at this 
problem by using the behaviour of the underlying hydrody- 
namical wind models - the sa me, which were used for the 
mass-loss parameterisation in iPhillips et all (|20Qgh - with 
respect to the outflow velocity. We are checking the assump- 
tion of constant outflow velocity in our models by checking 
the average velocities of the dust-driven wind models against 
dependence on the basic stellar quantities. 



2 DESCRIPTION OF THE WIND MODELS 

For this investigation we use a set of hydrodynamical wind 
models that include formation and growth of dust grains 
which, by radiation pressure, drive massive outflows. These 
are characteristic of the ultimate stages of stellar evolution 
on the AGB. 

The wind model s have been obta i ned w it h the computer 
code developed by iFleischer et al.l |l992l '), IWinters et al.l 
(2000), and references therein. They result from the selfcon- 



sistent solution of the non-linearly coupled system of equa- 
tions describing the hydrodynamical and thermodynami- 
cal structure of a spherically symmetric stellar atmosphere 
(with a pulsating photosphere as an inner boundary con- 
dition, providing an initial mechanical energy input), its 
chemical composition, as well as the nucleation, growth, and 
evaporation of dust grains, and the radiation of the central 
star. 

The hydrodynamical wind structure (mass density p 
and outflow velocity v) follows from the equation of continu- 
ity and the equation of motion which includes the radiation 
pressure on dust grains. The law of energy conservation and 
radiative transfer determine the temperature structure. 

A carbon-rich chemistry is assumed, where oxygen is 
completely locked in the CO molecule. The molecular com- 
position is calculated under the assumption of chemical equi- 
librium. The formation, growth, and evaporation of carbon 
grains is calculated accord i ng to the moment metho d deve l- 
oped bylGail fc Sedlmavrl l|l988h and lGaueer et al.1 (|l990h . 
ISchirrmacher et~al. ( 20031 ') give a more detailed summary of 
the physical assumptions. 

From a set of these models we derived a mass-loss for- 
mula for stars with solar ele ment abundances, w hich have 
reached the tip of the AGB |Wachter et al.ll2002h . For that 
set of wind models, the velocity amplitude of the piston, 
used to simulate stellar pulsation, was set to a value of 
5 km s _1 . That value was found to be in accordance with 
a sample of observed C-star lightcurves. The dependence of 
the mass-loss rates on the pulsation period was implicitly 
accounted for by the use of an observationally determined 
period-luminosity relation for Mira stars, the most appropri- 
ate class of objects on the tip of the AGB to compare with 
the hydrodynamical models. In this way we were able to 
represent the mass-loss rate by a simple relation dependent 
on the stellar parameters mass M, luminosity L, and effec- 
tive temperature T e ff, only. In this work, we follow a similar 
strategy in oder to find such a simple relation to represent 
the average wind velocity. 

The model set used for the current investigation is listed 
in Table [1] given are the input parameters stellar mass M, 
effective temperature T e g, luminosity L, carbon-to-oxygen 
ratio C/O, and piston period P. The resulting quantities av- 
eraged over typically 20 periods are mass- loss rate (M) , out- 
flow velocity {i> CX p), dust-to-gas ratio ( ^ dust ) , and radiative- 
to-gravitational acceleration ratio ( ^ rad ). 

The carbon-to-oxygen ratio of these models is 1.3, which 
is in good agree ment with both, recent stellar evolution 
calculations (see l| Weiss fc Ferguson! 120091 ) . Fig. 1) and the 
upper end of the observe d C/O range (see for instance 
iBergeat fc Chevallierl2005l ). According to both accounts, the 
C/O ratio rises with the final phases of stellar evolution, 
from just above 1 when entering the carbon star phase, to 
an extreme of around 1.3 on the very tip of the AGB, where 
we apply our dust-driven wind models. 



3 BEHAVIOUR OF THE WIND OUTFLOW 
VELOCITY 

Following simple physical arguments, the dynamics of op- 
tically thick winds should primarily depend on L/M, 
since this ratio is proportional to the above-mentioned 
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Figure 1. Outflow velocities for our selection of wind models with carbon-to-oxygon ratio C/O = 1.3; the solid line shows relation (JTJ. 



radiative-to-gravitational acceleration ratio { ° rad ). In fact, 
|Schroder et al.lll999h find that truly dust-driven outflows 
require an Eddington-like lower limit in the form of a critical 
L/M, above which the mean ratio of radiative to gravita- 
tional acceleration becomes (a) > 1. 

Motivated by this example of the relevance of L/M, 
Fig.[T]depicts the time-averaged outflow velocities as a func- 
tion of L/M, using the models detailed in Table[T] which are 
restricted to cases of dust-driven mass loss (characterized by 
(a) > 1, see above). Hence, Fig. [1] supports that there is a 
relation between the average wind velocity and the physical 
parameter of L/M. 

Based on this approach, we used the non- 
linear least square fit in the R-base package 
l|R Development Core Team! I2008T ) to derive the de- 
pendence of the outflow velocity on L/M in terms of a 
power function (L and M in solar units, « oxp in km s _1 ): 



/ r \ 0.57(±0.03) 

« oxp = 0.05(±0.02) (-) (1) 

This relation is depicted in Fig. [T] as a solid line. The 
relatively small exponent reflects the modest variation of 
the outflow velocity over the whole range of wind models 
depicted in Fig. [1] 



3.1 Velocity changes induced by stellar evolution 

In the next step, we model the changes of the mean wind 
velocity, which are induced by stellar evolution, driven by 
the related changes in stellar luminosity L and mass M on 
the tip of the AGB. In this context, a most critical case is 
the immediate aftermath of a thermal pulse (TP), where 
velocity changes could be suspected to occur on timescales 
rapid enough to compete with the dynamical timescale of 
the circumstellar envelope (several thousand years). For the 
purpose of studying this problem in more detail, we com- 
bine the above velocity-relation fl} with stellar evolution 
models, which consider the respective mass-loss of our dust- 
driven wind models at each time-step and yield L and M, 

accordingly. 

As in earlier work (see IPhillips et al 1 l2009l . and refer- 
ences therein) we use the well-calib rated, fast stella r evo- 
lution code originally developed by lEggletonl fl973). The 
parameterized mass- loss description applied in this code 
ijWachter et al.ll2002l ) is based on the same set of wind mod- 
els as the current investigation. As an example, Fig.[2]shows 
the last 80 000 years of mass-loss history of a 2.25 Mq star 
with solar element composition where we can see the super- 
wind phase on the tip of the AGB with a duration of ~20 000 
years. 

For the same model and timespan, we show the varia- 
tion of the outflow velocity in Fig. [3] The mean wind ve- 
locity here is seen to grow very slowly and steadily, with 
exceptional interruptions caused by thermal pulses. In these 
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Table 1. Set of Berlin wind models used in this investigations. 



The piston velocity amplitude is fixed to Av 
for details. 



5 km s , see text 
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Figure 2. Mass-loss history of a star with solar element abun- 
dances of initial mass M\ = 2.25 Mq. 




Figure 3. History of outflow velocity according to relation JTJ of 
same star as in Fig. [2] 



cases, the most prominent change is a temporary decrease, 
which coincide with the steep decrease in stellar luminosity, 
related to the temporary extinction of the hydrogen burning 
shell after a helium shell flash. The energy generated in the 
latter takes more time to reach the surface. It generates the 
slow recovery and further rise of luminosity, mass-loss rate 
and wind velocity until the next thermal pulse occurs. In 
any case, the relative variation of the wind velocity is very 
small compared to that of the large simultaneous changes in 
the mass loss rate (see Fig. And, important for hydro- 
dynamic considerations: No large, fast increase of the wind 
velocity, which would lead to a compression of the outflow 
ahead as suggested by colliding wind scenarios, is seen in 
Fig. O 

When considering the superwind phase as a whole, the 
wind velocity increases mostly steadily from 8 to 14 km s _1 
over a timespan of about 50 000 years. This is a lot longer 
than the dynamic age of any present-day outer cool PN enve- 
lope, which would contain evidence for only the past 10 000 
years or so. Over that latter order of timespan, only a very 
modest, systematic increase by about 1 km s _1 (about 10%) 
is expected according to Fig. [3] 
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Figure 4. Zoom on last thermal pulse of Fig. \3\ 




Figure 6. Comparison with wind models by iMattsson et al.l 
j201Ch with C-O = 8.5 (corresponds to C/O ~ 1.7); circles and 
triangles down: At) = 4 km s — 1 , triangles up and diamonds: 
Av = 6 km s" 1 , M = 0.75 and 1 Mq, respectively. 
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Figure 5. Temporal change over 40 periods of the outflow ve- 
locity in a hydrodynamical wind model - input parameters: 
M = 1 Mq, T off = 2600 K, L = 10000 L , P = 640 d, C/0=1.3, 
Ad = 5 km s —1 . 



3.2 The perspective: intrinsic outflow variation 

As the most rapid and largest velocity change, we identify 
the about 200 years after a thermal pulse (TP), when the lu- 
minosity goes through a rapid dip. The outflow velocity cor- 
responds here with a reduction of about 30%, on a timescale 
of only 50 years - see Fig. 21 which shows the final TP of 
Fig. |3] over a much smaller timespan. 

This is, in expansion timescale, still a modest and, yet, 
not very fast velocity change. It does not exceed much the 
intrinsic, erratic fluctuations of the outflow of the dynamic 
wind models, of which Fig.[5]shows a typical example. These 
frequently reach or exceed 10%, on a timescale of only 10 
years. We should therefore not expect to find any observable 
evidence for the TP-related changes of the wind velocity in 
the density structure of cool PN envelopes. Rather should 
the much larger relative changes of the mass-loss rate be 
held accountable for any distinguishable radial structure. 



4 DISCUSSION 

How do our findings compare with results from other dust- 
driven wind models? Recently, IMattsson et al.l ((2010) pub- 
lished a grid of wind models with solar element abundances. 
The major difference of their models is the inclusion of 
non-grey radiative transfer. For comparison, Fig. [6] shows 
the o utflow velocities of those models of iMattsson et al.l 
(2010), which have a free carbon abundance C-O = log[(nc — 
no)/nu] + 12 — 8.5 equivalent to a carbon-to-oxygen ratio 
of C/O ~ 1.7. This value is the best option to compare to 
our choice of C/O, since only a few models in their grid 
with C/O ~ 1.35, which is closer to our value, result in an 
outflow. Their L/M ratio does not extend to high values 
because in their grid models with luminosities higher than 
10 000 Lq have only been calculated for masses larger than 
1 Mq. 

Furthermore, the models of lMattsson et all (|2010l ) show 
a wider spread in their outflow velocities - possibly due 
the influence of effective temperature. As for the mechan- 
ical energy input, we have set the piston velocity ampli- 
tude to 5 km s - 1 . He nce, for Fig. [6] we used the models of 
IMattsson et al.l (|20 10f > with the closest choices in this re- 
spect, 4 and 6 km s _1 . The largest outflow velocity values 
seen in the figure correspond to the higher piston velocity. 

Even though their model set seems to suggest a steeper 
exponent - as would our models restricted to the lower L/M 
range - our velocity relation overall is consistent with their 
grid, especially taking int o account the wider spr ead of out- 
flow velocity values in thefMattsso n et al.l l|2010l ) models. 



5 CONCLUSIONS 

From the analysis of a consistent set of dust-driven wind 
models, we find that the velocity changes of the superwind 
can be parameterised by a simple dependence on L/M. In 
absolute terms, the related velocity changes are of a small 
nature, of the order of a fraction of the typical wind velocity. 
This is consistent with the very modest dependence of the 
velocity on L and M found by lHabing et al.l (|l994), and 
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both these findings are contrasted by the huge changes of 
the mass-loss rate (by several orders of magnitude) . 

Looking at the final phase of the superwind, systematic 
velocity changes are not exceeding the inherent short-term 
fluctuations of the dust-driven wind. And when a TP falls 
into this phase, there should only be a 30% velocity dip 
on a timescale of 50 years, which should still be difficult 
to extract from the density profile in view of the intrinsic 
erratic fluctuations of a dust-driven wind (over f0% on a 
timescale of 10 years). 

Hence, we may conclude that the huge changes of the 
mass-loss rate in the history of a superwind, especially in 
the aftermath of a TP, are by far the dominant factor in 
shaping a present-day cool PN envelope, while changes in 
the velocity are hardly exceeding the natural fluctuation of 
a dust-driven wind. This find gives some support to a inter- 
pretation of cool PN envelopes in terms of a simple, constant 
outflow history, where the density profile mostly shows the 
progressively (in radial direction, that is with outflow age) 
diluted superwind and its mass-loss rate evolution. 
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